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We report the results of a new accurate evaluation of light nuclei yields in primordial nucleosynthesis. The 
relic densities of 4 He, D and 7 Li have been numerically obtained via a new updated version of the standard BBN 
code. 



1. Introduction 

Big Bang Nucleosynthesis (BBN) is one of the 
most powerful tools to study fundamental inter- 
actions since light nuclei abundances crucially de- 
pend on many elementary particle properties, like 
the number of effective neutrino degrees of free- 
dom, N v . At the moment, however, recent data 
on He mass fraction, Y4, and Deuterium (D) and 
7 Li abundances, Y2 = D/H and Yj = 7 Li/H, 
produced during BBN are controversial, since 
there are different sets of results, two of them mu- 
tually incompatible: Y} 1) = 0.234 ±0.002 ±0.005 



and Y, 



(h) 



0.243 ± 0.003 InJ, Yr 
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in 



0.3)10" 5 and F 2 ( ' = (1.9 ± 0.4)10~ 4 
(1.6 ± 0.36)10- 10 and Y 7 (h) = (1.73 ± 0.21)10" 
H (for a brief summary of the experimental sit- 
uation on primordial abundances see Ref. H). 

In spite of this discrepancies, which could be of 
systematic origin, the above results for A He data 
indicate that one is reaching a precision of the 
order of percent, requiring a similar level for the 
uncertainties on the theoretical predictions. Be- 
sides the corrections to the proton/neutron con- 
version rates, which fix at the freeze out temper- 
ature ~ 1 MeV the neutron to proton density 
ratio, the other main source of theoretical uncer- 
tainty comes from nuclear rates relevant for nuclei 
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formation. In some cases, these rates are known 
to well describe the data in a temperature interval 
which only partially overlaps the one relevant for 
BBN, 0.01 MeV < T < 10 MeV. Recent studies 
||, however, show that, in particular for the 4 He 
mass fraction, the effect is at most as large as the 
one due to the uncertainty on neutron lifetime 
r n , and smaller than 1%. Therefore it is theoreti- 
cally justified to look for all sources of theoretical 
uncertainty up to this level of precision. 

In a previous paper Q we performed a thor- 
oughly analysis of all corrections (electromagnetic 
radiative corrections, finite nucleon mass correc- 
tions, plasma effects) to the rates of the processes 
converting n <-> p, i.e. v e n «-> e~ p, T> e p «-> e + n 
and n <-> e~ V e p. Here, we report on a follow- 
ing work Q), where we included the above men- 
tioned corrections in a new updated version of 
the standard BBN code ||. This new code was 
used for integrating the set of equation of BBN 
and obtaining the values of the primordial light 
nuclei yields. From the comparison of these pre- 
dictions with the experimental abundances it is 
possible to get informations on the effective num- 
ber of neutrinos and the final baryon to photon 
density ratio, 77. 

2. Corrections to Born rates 

As is well known, the key parameter in deter- 
mining the primordial 4 He mass fraction is the 
value of the neutron to proton density ratio at the 
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freeze-out temperature T ~ 1 MeV , since almost 
all residual neutrons are captured in £/e nuclei 
due to its large binding energy per nucleon. We 
shortly summarize the kind of corrections which 
are studied in detail in ||. 

The Born rates, obtained in the tree level V — 
A limit and with infinite nucleon mass have to 
be corrected to take into account basically three 
classes of relevant effects: 

i) order a radiative corrections. These effects 
have been extensively studied in literature 
and can be classified in outer factors, in- 
volving the nucleon as a whole, and inner 
ones, which instead depend on the details of 
nucleon internal structure. Actually, other 
small effects are expected at higher order in 
a, since the theoretical value of the neutron 
lifetime is compatible with the experimental 
one, 7%" = 886.7±f .9 s §], at 4-er level only. 
These additional contributions are usually 
taken into account by eliminating the cou- 
pling in front of the reaction rates in favour 

of T™. 

ii) All Born amplitudes should also be cor- 
rected for nucleon finite mass effects. They 
affect both the weak amplitudes, which 
should now include the contribution of nu- 
cleon weak magnetism, and the allowed 
phase space. Initial nucleons with finite 
mass will also have a thermal distribution 
in the comoving frame, producing a third 
kind of finite mass correction. 

iii) Since all reactions take place in a thermal 
bath of electron, positron, neutrinos, an- 
tineutrinos and photons, thermal-radiative 
corrections should be also included, which 
account for the electromagnetic interactions 
of the in/out particles with the surrounding 
medium. They can be evaluated in the real 
time formalism for finite temperature field 
theory. 

After solving the BBN set of equations, one can 
see that for all nuclides the pure radiative cor- 
rection provides the dominant contribution, while 
the finite nucleon mass effects and the thermal- 
radiative ones almost cancel each other. 



The total proton/neutron conversion rates were 
fitted, in the range 0.01 MeV < T < 10 MeV, to 
the following functional forms, 
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where z is the dimensionless inverse photon tem- 
perature, z = m e /T and the values of the param- 
eters can be found in Jjj . Note that Eq. (||) is valid 
only in the range 0.1 MeV < T < 10 MeV, be- 
cause oj p -> n ~ for T < 0.1 MeV. The fits have 
been obtained requiring that the fitting functions 
differ by less than 0.1% from the numerical values 
in the considered range. 

3. Numerical code 

3.1. The Equations of BBN 

Denoting with R the universe scale factor, ng 
the baryonic density, cf) e = HejT the electron 
chemical potential, and Xi the nuclide number 
densities, X{ = ni/ns, the BBN set of equations 
is a system of coupled differential equations in the 
previous unknown functions of time. By expand- 
ing the equations with respect to <p e and changing 
the evolution variable to z, after a little algebra 
one is left with the following N nuc + 1 equations, 



G(z,h,Xj) 



r, 



(3) 
(4) 



dh r , t , „, j 3/i 

— = 1 -H(z,h,Xj) G(z,h,Xj} 

dXj 
dz 

where we have introduced the dimensionless 
baryon density, h = ub/T 3 , Hubble parameter, 
H = H/m e , and nuclear rates T,; = Ti/m e . The 
function G in Eq.s (0) and M) is 



G(z,h,X 3 ) 



xp 



^2(Pa +Pa) 



H + h 
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Table 1 



The predictions on light element abundances obtained with the numerical code for 77 = 5 • 10 
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0.3638-10" 4 


0.1175-10" 4 


0.2446 


0.2814-10" y 




0.3727-10" 4 


0.1184-10" 4 


0.2550 


0.2873-10" 9 
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In the previous equation zn = m e (MeV)/2.3 
is the inverse neutrino decoupling temperature, 
a = e,7, M u = M u /m e and AMj = AM,/m e = 
(Mi — AiM M )/m e are the dimensionless atomic 
mass unit and mass excess. We have neglected, in 
the original system, terms containing the deriva- 
tives of chemical potential. The expression of 
the dimensionless energy densities and pressures, 
p Q = p a /T 4 and p a = p a jT 4 , contained in G, 
were evaluated taking also into account the 7 and 
electromagnetic mass renormalization, and fit- 
ted as functions of z for their inclusion in the BBN 
code (see Appendix A of (7]]). The previous effect, 
changing the 7 and e^- equations of state, modi- 
fies the T v /T ratio too. However, the difference 
between the neutrino temperature evaluated with 
the correct rcnormalized masses and the one ob- 



tained with approximated expressions, m\ 







and 



aT 2 /m e , results to be smaller than 



0.01%, a correction that can be neglected at the 
level of precision we are interested in. 

The initial conditions for Eq.s (0) and (|4|) are 
given by 
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in terms of the final baryon to photon density 
ratio, 77, and 
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which represents the condition of nuclear statisti- 
cal equilibrium for an arbitrary i-th nuclide, with 
gi internal degrees of freedom, Zi and Ai charge 



and atomic number, and Bi binding energy. This 
condition is satisfied with high accuracy at the 
initial temperature Tj n = 10 MeV. 

3.2. Numerical Method 

The numerical problem of solving the set of 
equations (||) and (|]) is stiff, because the r.h.s. of 
(0) results to be a small difference of large num- 
bers. While in the standard code || the implicit 
differentiating method (backward Euler scheme) 
fic| l for writing the r.h.s. of ((J) and a Runge- 
Kutta solver are used, we choose a method be- 
longing to the class of Backward Differentiation 
Formulas (BDFs) (l(J, implemented by a NAG 
routine. 

The new code includes all the 88 reactions be- 
tween the 26 nuclides present in the standard 
code with the same nuclear rate data, collected 
and updated in jnj. However, in order to reduce 
the computation time we used a reduced network, 
made of 25 reactions between the first 9 nuclides 
(see Table 1 of Ref. ) , verifying that this affects 
the abundances of light nuclei for no more than 
0.01%. 

4. Results and conclusions 



We report in Table Ft] the predictions for Y2 — 



M 6 Xr 



and Y7 



x 2 



responding to the complete n p rates, uJ r p i, 
and to the Born approximation, LOg. This last 
quantity denotes the pure Born predictions for 
n <-> p rates without any constant rescaling of 
coupling to account for the experimental value 
of neutron lifetime. These values have been ob- 
tained for N v — 3 and r\ — 5T0 -10 . The net effect 
of the corrections is to allow a smaller number of 
neutrons to survive till the onset of nucleosynthe- 
sis. This ends up in a smaller fraction of elements 
fixing neutrons with respect to the pure hydrogen. 
In Fig. the predictions on Y4 are shown for 
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N u = 2,3,4 and for a 1 a variation of t^ x . The 
three solid lines are, from larger Y4 to lower val- 
ues, the predictions corresponding to N v = 3 
and = 888.6 s, 886.7 s, 884.8 s, respec- 
tively. Analogously, the dashed lines correspond 
to N v = 4 and the dotted ones to N v = 2. The 
experimental estimates, as horizontal bands, are 
also reported. In Fig.s || and || the D and 7 Li 
abundances are reported with the same notation. 
Note that, due to the negligible variation of Y2 
and Y7 on small r n changes, no splitting of pre- 
dictions for 1 a variation of t^ x is present. 

By fitting, up to one percent accuracy, Y2, ¥3, 
Y4 and Yj as a function of x — log 10 (lO 10 77), 
N v and r„, the following expressions have been 
obtained: 



Figure 2. The quantity Y2 versus rj is reported. 
The horizontal dashed and dotted bands are the 
experimental values. 



io 3 -y 2 = 



<H x i + a 5 (N v - 3) 
+a 7 x 2 } , 



exp {— x 



(8) 



Figure 4. The likelihood distributions for the light element yields Y2, I4, Y7 are shown as functions of 
N v and log 10 10 10 ?7, in arbitrary units. From left to right and from top to bottom the following cases are 
considered: a) high D, low 4 He; b) high D, high 4 He; c) low D, low 4 He; d) low D, high 4 He. The 
normalizations of c) and d) are 25 and 100 times, respectively, the one of a) and b). 
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io b -y = 
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a, x i + a 5 (N u - 3) 



exp{-a 6 x} ,(9) 



10-y 4 = ^ a, x l + a 6 (r - r ex ) + a 7 (N v - 3) 



i=0 



io y -y 7 = 



+a 8 x (t - r ex ) + a 9 x (iV„ - 3), 

3 

^ Oi x* + a 4 (iV„ - 3) + a 5 x (iV v - 3) 



(10) 



x exp {— etg x + a 7 x 2 + a% x 3 + ag x 4 } , (11) 

where the values of the fit coefficients are reported 
in Table 3 of §. 

In Fig. [I] we plot the product of gaussian dis- 
tribution for D, A He and 7 Li, centered around 
the measured values and with their correspond- 
ing experimental errors, 



L{N v ,x) 



exp 



x exp 



x exp 



(Y 2 (N U) x)- 




2a 2 




(Yi(N v ,x)- 


yexy 






(Y 7 (N u ,x)- 
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2a 2 7 



(12) 



The experimental values used in the previous 
equation correspond to the four combinations of 
experimental results: a) y£ , Y^; b) Y% , 3^ ; 

c) y a w, yw ; d) y a w, yf ). 

The figure shows that the high value of D is 
preferred (plots a) and b)). In both cases the 
distributions are centered in the range x 6 0.2 

0. 4, but at N v ~ 3 for low i He and N„ - 3.5 
for high 4 He. For low D the compatibility with 
experimental data is worse (note that in c) and 

d) cases the distributions have been multiplied by 
a factor of 25 and 100 respectively) and centered 
in the range x G 0.6 -j- 0.8, and at N v ~ 2 for low 
4 #e and iV^ ~ 3 for high A He. 
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